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INTRODUCTION 



This subroutine set provides an IMP-16 microprocessor with double-precision and floating-point capability. 
Double-precision and floating-point computation is very useful when additional precision is needed or the 
range of numbers is expected to be large. In the following pages, individual subroutines are described alon§ 
with necessary user instructions. Subroutines included in the set are listed in the table below. 



ARITHMETIC SUBROUTINE SET 



SUBROUTINE 


MNEMONIC 


ENTRY POINT 


Single Precision Multiply 


MULT 


FCOl 


Single Precision Divide 


DIV 


FCIA 


Double Precision Multiply 


DPMUL 


FCCF 


Positive Double Precision Divide 


DPDIV 


FDE3 


Double Precision Square 


DPSQUARE 


FCCD 


Double Precision Complement 


DPCOMP 


FDBA 


Double Precision Shift 


DPSH 


FDDB 


Double Precision Shift Right 


DPSHR 


FDCE 


Double Precision Shift Left 


DPSHL 


FDDB 


Quadrant Tests 


QUAD 


FDAO 


Sine 


SIN 


FCFA 


Cosine 


COS 


FCFC 


Arctangent 


ARCTAN 


FD37 


Floating Point Add 


FPADD 


FC30 


Floating Point Multiply 


FPMUL 


FC5D 


Floating Point Divide 


FPDIV 


FC6C 


Floating Point Complement 


FPCOMP 


FCB9 


Check Zero Exponent 


CZERO 


FC63 


Extract Exponent to Stack 


EXTEXP 


FC7C 


Add Exponent from Stack 


ADDEXP 


FCBD 


Left Normalize 


LFNOR 


FC97 


Double Left Normalize 


DLNORM 


FC8A 


Fraction to Floating Point 


FLOAT 


FC85 


Floating Point to Fraction 


SFO 


FCC7 



GENERAL CHARACTERISTICS 



The arithmetic subroutine set uses 512 memory words located at addresses FCOO through FDFF (64512 
through 65023, decimal). All subroutines are contained wholly within the set; they require no external 
program code. The IMP-16 extended instruction set is required. Operands and results, in most cases, 
use the general registers. The only memory locations written into by the set are in the range OOEO 
through OOEF. All constants are contained within the subroutine set. 



EXECUTION TIMES 



Execution times shown on each of the subroutine data sheets aTe expressed in terms of three variables, where 



R = Number of main memory read cycles 
W = Number of main memory write cycles 
N = Number of microprogram cycles 



These execution time expressions can thus be used to calculate subroutine execution time for any IMP-16 
implementation. A few subroutine execution times are also dependent on user -specified variables (such as 
shift counts in shift instructions). These are defined on the individual data sheets. All execution times are 
based on worst -case conditions. 

NUMBER REPRESENTATION 

Numbers in a computer program often represent real world parameters. In this subroutine set, four different 
modes of numerical representation are used: integer, fractional, double-precision, and floating-point. All 
four modes have important applications in IMP-16 systems. 

Integer notation is the most common and obvious mode of numerical representation. It is simply an integer 
count of something. In the IMP-16, numbers up to 2^^ may be represented because of the 16-bit word length. 
To represent both positive and negative integers, a twos -complement system is used. This gives an integer 
range of 

-2^^ < integer < 2^^ 

With twos -complement notation, the most significant bit, bit 15, can be thought of as a sign bit, and the binary 
point can be thought of as being to the right of bit 0. 

FRACTIONAL NOTATION 

Many real -world parameters are continuous rather than discrete. Such values can never be indicated with per- 
fect precision because of the limited word length of any computer. Take the airspeed of an aircraft as an 
example. It may have a potential range from zero to 1000 miles per hour. We would rarely expect its speed 

to be an exact round number. At some instant the speed may be 556. 234 If we want to express this 

speed as precisely as possible given a 16-bit word length, an integer notation of 556 miles per hour would not 
be the best choice. A more convenient technique is to use fractional notation. The parameter, in this case 

speed, can be scaled so that 1000 miles per hour is equivalent to 1. 0000 This woiild allow the full 16-bit 

capacity to be utilized. After scaling, the range of numerical expression of a fractional number is 

- 1.0 < N < 1.0 

With fractional notation, the binary point can be thought of as being to the right of bit 15. Thus, the largest 
positive binary number that can be represented is 0. 11111. . . or almost one. Precision is about 3.1* 10"^. 

DOUBLE PRECISION 

When precision greater than 16-bits is required, double-precision numbers can be used. Double-precision 
numbers give a precision of one part in 2 or about 4. 6* lO"-*^^. Double-precision numbers are an extended 
fractional notation. Note that two IMP-16 words are required. 

FLOATING POINT 

Floating-point representation is a technique to express numbers in a form similar to scientific notation, with 
a fractional mantissa and an exponent. The value of any floating-point number so represented is equivalent to 
mantissa * (2**exponent). The mantissa consists of 24 bits and gives a precision of one part in 2^^. The 
exponent has 8 bits and allows exponents in the range 

7 7 

-2 < exponent < 2 



The great advantage of floating-point numbers is that the scaling factor of a fractional number becomes a 
visible attribute which can be manipulated. Thus, the scaling factor can be altered to maintain precision 
or facilitate computations. 



NUMBER REPRESENTATION EXAMPLES 
li 



INTEGER 



15 



-?^ 







BINARY POINT 



EXAMPLES: 

5 = 000516 
-5 = FFFB 



16 



RANGE: -32768< N < 32768 
MAXIMUM PRECISION: 1 
SMALLEST VALUE: 1 



SPEED = 556. 234 = 556 = 022C, ^ 
10 16 



• FRACTIONAL 



-Ih 



15 . 14 



\ 



BINARY POINT 

RANGE: -1 < N < 1 

MAXIMUM PRECISION: 1/2^^ = 3.1 * 10 "^ 

SMALLEST VALUE: 2"^^ 



EXAMPLES: 
1/2 = 4000j^g 

-l/2=C000i6 
-1 = SOOOie 

556.234... 
SPEED = 1000.000 



= 0. 556234^Q = 4732 



^16 



DOUBLE PRECISION (FRACTIONAL) 



31. 30 , , 



RANGE: -1 < N < 1 

MAXIMUM PRECISION: l/2^^ = 4. 6 * 10"^^ 

SMALLEST VALUE: 2"^^ 



16 



15 



EXAMPLES: 
1/2 = 4000 OOOO^^g 
-1/2 =C000 0000 



-1 = 8000 0000 



16 



16 



SPEED = 4732 ACFBj^g 



FLOATING POINT 



RANGE: -2^^'^< N< 2^^^ 
PRECISION: 2~^^ (g^^P^^^N"^) 
SMALLEST MAGNITUDE: 2~^^^ 
LARGEST MAGNITUDE: (1-2"^^) 2^^"^ 



-&- 



23. 22 



^ 



MANTISSA (M) 




EXPONENT (E) 



N = M * 2^ 

EXAMPLES: 
1/2 = 4000 00i6 (00i6) 
1/2 = 2000 00i6 (Olie) 
1/2 = 1000 00i6 (02i6) 
1/4 = 2000 00i6 (OOie) 
1/4 = 4000 00i6 (FFi6) 



SPEED = 556. 234 = 4587 TA-^q (OA^g) 



ANGULAR REPRESENTATION 



Aagles are expressed as a double -precision fraction of tt radians. 



0000 



EOGO 



COOO 



16 



AOOO 




4000 



16 



NOTATION 

In the descriptions which follow, RO, Rl, R2, and R3 represent the contents of the four accumulators = AGO, ACl, 
AG2, and ACS. 



Normalization of a floating point number means that the mantissa is shifted left until the most significant bit is a 1, 
and the exponent adjusted accordingly. 



SUBROUTINE: SINGLE PRECISION MULTIPLY 



LABEL: MULT 



DESCRIPTION: MULTIPLIES TWO SINGLE -PRECISION, SIGNED, FRACTIONAL 16-BIT NUMBERS TO 
GIVE A SIGN PLUS A 30-BIT RESULT 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FCOl 










R0=X 


R2=Y 








Rq, r^ = 


= X * Y R2 


- ALTERED 


R3 -- 


= UNCHANGED 



30R + 243N + W 



TEMPORARY LOCATIONS USED: 
SUBROUTINES CALLED: 
STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



OOEO 

DSHL, DPCOMP 

2 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

2. FRACTIONAL DIFFERS FROM AN INTEGER MULTIPLY IN THAT THE RESULT OF A 
FRACTIONAL MULTIPLY IS SHIFTED LEFT 1 BIT. 



SUBROUTINE: SINGLE PRECISION DIVIDE 



LABEL: DIV 



DESCRIPTION: DIVIDES A SIGNED, FRACTIONAL 16-BIT NUMBER INTO A SIGNED, FRACTIONAL 
32-BIT NUMBER, PROVIDING A QUOTIENT AND AN INTEGER REMAINDER 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FCIA 
Rn, R 



= DIVIDEND 



Ro = DIVISOR 



Rq = QUOTIENT R^ = REMAINDER 
49R + W + 343N 



R2 = ALTERED R3 



UNCHANGED 



TEMPORARY LOCATIONS USED: 
SUBROUTINES CALLED: 
STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



OOEO 

DSHR, DPCOMP 

2 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
2. OVERFLOW FLAG IS SET IF MSB OF DIVIDEND > DIVISOR. 



SUBROUTINE: DOUBLE PRECISION MULTIPLY 



LABEL: DPMUL 



DESCRIPTION: MULTIPLIES DOUBLE -PRECISION OPERANDS TO GIVE A DOUBLE -PRECISION PRO- 
DUCT 



ENTRY POINT: FCCF 

ENTRY CONDITIONS: Rq, R^ - X Rg, R3 = Y 

EXIT CONDITIONS: Rq, R^ = X * Y Rg, R3 ALTERED 

EXECUTION TIME: 200R + 8W + 1125N 



TEMPORARY LOCATIONS USED: OOEO, OOEl, OOEE, OOEF 

SUBROUTINES CALLED: QUAD, MULT, DPCOMP, DSHL 

STACK WORDS NEEDED: 4 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: POSITIVE DOUBLE PRECISION DIVIDE 



LABEL: DPDIV 



DESCRIPTION: UNSIGNED FRACTIONAL DIVIDE OF Ro, R3 INTO Ra, R,, BOTH DOUBLE -PRECISION 
POSITIVE FRACTIONS 

ENTRY POINT: FDE3 

ENTRY CONDITIONS: Rq, R^ = POSITIVE DIVIDEND Rg, R3 = POSITIVE DIVISOR 

EXIT CONDITIONS: Rq, R^ = QUOTIENT R2, R3 - REMAINDER SELFF = 

EXECUTION TIME: 510R + 97W = 2578N 

TEMPORARY LOCATIONS USED: OOEO, OOEC-OOEF 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: DOUBLE PRECISION SQUARE 



LABEL: DPSQUARE 



DESCRIPTION: SQUARES THE DOUBLE -PRECISION NUMBER IN Rq, R-,^ 



ENTRY POINT: FCCD 

ENTRY CONDITIONS: Rq, R;,^ = X 

EXIT CONDITIONS: Rq, R;^ = X^ Rg, R3 ALTERED 

EXECUTION TIME: 202R + 8W + 1137N 



TEMPORARY LOCATIONS USED: OOEO, OOEl, OOEE, OOEF 

SUBROUTINES CALLED: DPMULT, QUAD, MULT, DPCOMP, DSHL 

STACK WORDS NEEDED: 4 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSED ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: DOUBLE PRECISION COMPLEMENT 



LABEL: DPCOMP 



DESCRIPTION: COMPUTES TWOS COMPLEMENT OF A DOUBLE -PRECISION NUMBER 



ENTRY POINT: FDBA 

ENTRY CONDITIONS: Rq, R^^ - X 

EXIT CONDITIONS: Rq, R-^ = TWOS COMPLEMENT OF X Rg, R3 PRESERVED SELFF = 

EXECUTION TIME: 16R + 53N 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: DOUBLE PRECISION ARITHMETIC SfflFT 



LABEL: DPSH 



DESCRIPTION: SHIFTS A DOUBLE -PRECISION NUMBER TO RIGHT OR LEFT AS SPECIFIED BY Rg. 
ABSOLUTE VALUE OF Rg INDICATES NUMBER OF BITS TO SHIFT. POSITIVE Rg 
INDICATES LEFT SHIFT; NEGATIVE Rg MEANS RIGHT SHIFT. 

ENTRY POINT: FDDB 

ENTRY CONDITIONS: Rq, Rj^ = X Rg = M (SHIFT COUNT) 

EXIT CONDITIONS: Rq, R-^ - SHIFTED X Rg = ALTERED R3 PRESERVED 

EXECUTION TIME: LEFT 5R + 20N + [SR + 32N] M 

RIGHT 7R + 26N + [4R + 21N] M 

TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

RIGHT SHIFT IS ARITHMETIC; THAT IS, SIGN BIT IS PRESERVED 



SUBROUTINE: DOUBLE PRECISION ARITHMETIC SHIFT RIGHT 



LABEL: DPSHR 



DESCRIPTION: SHIFTS DOUBLE -PRECISION NUMBER RIGHT M BITS, WHERE Rg = M 



ENTRY POINT: FDCE 

ENTRY CONDITIONS: Rq, R^ = X Rg = M (BITS TO SHIFT) 

EXIT CONDITIONS: Rq, Ri = RIGHT SHIFTED X Rg ALTERED R3 PRESERVED 

EXECUTION TIME: 5R + 20N + (4R + 21N) M 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

SIGN BIT IS PRESERVED 



SUBROUTINE: DOUBLE PRECISION SHIFT LEFT 



LABEL: DPSHL 



DESCRIPTION: SHIFTS DOUBLE -PRECISION NUMBER LEFT M BITS, WHERE Rg = M 



ENTRY POINT: FDDB 

ENTRY CONDITIONS: Rq, R;^ = X Rg = M (NUMBER OF BITS TO SHIFT) 

EXIT CONDITIONS: Rq, R^ = LEFT SHIFTED X R2 ALTERED R3 PRESERVED 

EXECUTION TIME: 5R + 20N + (8R + 32N) M 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: QUADRANT TESTS 



LABEL: QUAD 





^1 


^2 





+ 


+ 


1 


+ 




2 


- 


- 


3 


- 


+ 



DESCRIPTION: GIVEN TWO SIGNED DOUBLE -PRECISION NUMBERS CI AND C2, THE INITIAL SIGN OF 
EACH IS INDICATED BY A FLAG IN MEMORY AND THE ABSOLUTE VALUE OF THE 
NUMBERS IS GIVEN. 



ENTRY POINT: FDAO 

ENTRY CONDITIONS: Rq, R^ = C^ R2' ^3 " ^2 

EXIT CONDITIONS: Rq, R-^ = ABS (C^) Rg, R3 = ABS (C2) OOEl = 



EXECUTION TIME: 47R + W + 175N 



TEMPORARY LOCATIONS USED: OOEl 

SUBROUTINES CALLED: DPCOMP 

STACK WORDS NEEDED: 1 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



WHERE ZERO IS + 



SUBROUTINE: SINE 



LABEL: SIN 



DESCRIPTION: CALCULATES SINE OF DOUBLE -PRECISION ANGLE IN Rq, Rj 
CALCULATION IS MADE: SIN(X) = COS (X - •^/2) 
ACCURACY IS 0. 00000 0002y]^Qj 



ENTRY POINT: FCFA 

ENTRY CONDITIONS: R , R^ = X 

EXIT CONDITIONS: Rq, Rj^ = SIN(X) Rg, Rg ALTERED 

EXECUTION TIME: 1697R + 82W + 9300 N 



TEMPORARY LOCATIONS USED: OOEO-OOEB, OOEE, OOEF 

SUBROUTINES CALLED: DPCOMP, DPMULT, PEXPN 

STACK WORDS NEEDED: 6 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: COSINE 



LABEL: COS 



DESCRIPTION: CALCULATES COSINE OF DOUBLE -PRECISION ANGLE IN R^, Rj 
POLYNOMIAL APPROXIMATION USED IS: 

COS(X) = 1 + X^ (Ag + X^ (A4 + X^ (Ae + X2 (Ag + X^ (A^^q))))) 
WHERE: Ag = -0.49999 99963 Ag = 0. 00002 47609 

A4 = 0. 04166 66418 A^q = "0. 00000 02605 

Ag = -0. 00138 88397 

FOR ANGLES IN THE FIRST QUADRANT 
ACCURACY IS 0. 00000 0002^3^0^ 



ENTRY POINT: FCFC 

ENTRY CONDITIONS: Rq, R^ = X 

EXIT CONDITIONS: Rq, R^ = COS(X) Rg, Rg ALTERED 

EXECUTION TIME: 1677R + 82W + 9285N 



TEMPORARY LOCATIONS USED: OOEO-OOEB, OOEE, OOEF 

SUBROUTINES CALLED: DPCOMP, DPMULT, PEXPN 

STACK WORDS NEEDED: 6 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



10 



SUBROUTINE: ARCTANGENT 



LABEL: ARC TAN 



DESCRIPTION: COMPUTES ARCTANGENT (Cl/C2), WHERE CI AND C2 ARE DOUBLE -PRECISION 
NUMBERS IN THE RANGE: 

-1 < C1,C2 < 1 

THE POLYNOMIAL APPROXIMATION USED WAS: 

ARCTAN (X) = X (1 + X^ (Ag + X^ (A4 + X^ (Ag + X^ (Ag + X^ (A^q + X^ (A^^g "^ ^^ 
(A^4 + x2 (A^g)))))))) WHERE 

A2 = -0. 33333 14528 A^q = -0. 07528 96400 

A4= 0.19993 55085 ^2^ 0.04290 96138 

Ag = -0. 14208 89944 ^lA = "0- 01^16 57367 

Ag = 0. 10656 26393 A-^g = 0. 00286 62257 

ACCURACY IS 0. 00000 002(io) 



ENTRY POINT: FD37 

ENTRY CONDITIONS: Rq, R^ = C;^ Rg, R3 = Cg 

EXIT CONDITIONS: Rq, R^ = ARCTAN (C^^/Cg) Rg. R3 ALTERED 

EXECUTION TIME: 2985R + 231W + 15892N 



TEMPORARY LOCATIONS USED: OOEO-OOEF 

SUBROUTINES CALLED: QUAD, DPCOMP, DPDIV, DPMUL, PEXPN, DSHR 

STACK WORDS NEEDED: 7 
(MAXIMUM DEPTH) 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

2. SINGULARITY CASE WHERE CI = C2 = YIELDS ARCTAN(Cl/C2) = 0. 



SUBROUTINE: FLOATING POINT ADD 



LABEL: FPADD 



DESCRIPTION: ADDS TWO FLOATING-POINT NUMBERS. 



ENTRY POINT: FC30 

ENTRY CONDITIONS: Rq, R^=X Rg, R3 = Y 

EXIT CONDITIONS: Rq, R^ = X + Y Rg, R3 ALTERED 

EXECUTION TIME: 147R + 7W + 613N 



TEMPORARY LOCATIONS USED: OOEO, OOEl, 00E2 

SUBROUTINES CALLED: EXTEXP, DSHR, CZERO, ADDEXP 

STACK WORDS NEEDED: 1 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: FLOATING POINT MULTIPLY 



LABEL: FPMUL 



DESCRIPTION: MULTIPLIES TWO FLOATING-POINT NUMBERS 



ENTRY POINT: FC5D 

ENTRY CONDITIONS: R^, R^ = X Rg, R3 = Y 

EXIT CONDITIONS: Rq, Rj^ = X * Y Rg, Rg ALTERED 

EXECUTION TIME: 1215R + 63W + 5077N 



TEMPORARY LOCATIONS USED: OOEO, OOEl, 00E5, 00E6, OOEE, OOEF 

SUBROUTINES CALLED: DLNORM, DPMUL, ADDEXP, CZERO 

STACK WORDS NEEDED: 4 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: POSITIVE FLOATING POINT DIVIDE 



LABEL: FPDIV 



DESCRIPTION: DIVIDES TWO POSITIVE FLOATING-POINT NUMBERS 



ENTRY POINT: FC6C 

ENTRY CONDITIONS: Rq, R^^ = X Rg, R3 = Y 

EXIT CONDITIONS: Rq, R^ = X/Y 

EXECUTION TIME: 1540R + 152W + 6584N 



TEMPORARY LOCATIONS USED: OOEO, OOEl, 00E2, 00E5, 00E6, OOEC, OOED, OOEE, OOEF 

SUBROUTINES CALLED: DLNORM, DSHR, DPDIV, ADDEXP, CZERO 

STACK WORDS NEEDED: 4 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: FLOATING POINT COMPLEMENT 



LABEL: FPCOMP 



DESCRIPTION: COMPLEMENTS A FLOATING-POINT NUMBER IN R^, R^^ 



ENTRY POINT: FCB9 

ENTRY CONDITIONS: Rq = Xi R-^ = Xg, EX 

EXIT CONDITIONS: Ra = ^ Ri = ^, EX R2, R3 PRESERVED 

EXECUTION TIME: 50R + 2W + 183N 



TEMPORARY LOCATIONS USED: OOEO 

SUBROUTINES CALLED: EXTEXP, DPCOMP, ADDEXP 

STACK WORDS NEEDED: 1 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: CHECK ZERO EXPONENT 



LABEL: CZERO 



DESCRIPTION: IF MANTISSA OF A FLOATING-POINT NUMBER IS ZERO, EXPONENT IS FORCED TO 
ZERO. 



ENTRY POINT: FC63 

ENTRY CONDITIONS: Rq, R^ = Xi, Xg (EX) 



EXIT CONDITIONS: IF Xi, Xg = Rq, R^ = 0,0(0) ^ PRESERVED 

IF X^, X2 7^ Rq, R2 = X^, X2 (EX) J 
EXECUTION TIME: 36R + 2W + 139N 



TEMPORARY LOCATIONS USED: OOEO 

SUBROUTINES CALLED: EXTEXP, ADDEXP 

STACK WORDS NEEDED: 1 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: EXTRACT EXPONENT TO STACK 



LABEL: EXTEXP 



DESCRIPTION: EXTRACTS EXPONENT FROM A FLOATING-POINT NUMBER, STORES EXPONENT ON 
STACK, AND CLEARS EXPONENT FIELD OF NUMBER. 



ENTRY POINT: PC7C 

ENTRY CONDITIONS: Rq, Rj^ = Xj^, Xg (EX) 

EXIT CONDITIONS: Rq, R^ = Xj^, Xg (0) STACK = EX 

EXECUTION TIME: 13R + W + 51N 



Rg, R3 PRESERVED 



TEMPORARY LOCATIONS USED: OOEO 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



COMMENTS: 



ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: ADD EXPONENT FROM STACK 



LABEL: ADDEXP 



DESCRIPTION: AN EXPONENT IN THE STACK IS PLACED IN A FLOATING-POINT WORD IN Rq, R^. 



Ro 




E 


1 




STACK 


Xl 




X2 


X3 




EX 


NEXT 


BEFORE 


' 


>A< 


^ 


^ 


Xl 




X2 


EX 




NEXTj 


NEXT2 


AFTER 















ENTRY POINT: FCBD 

ENTRY CONDITIONS: Rq, Rj^ = X^, Xg (X3) STACK = EX 

EXIT CONDITIONS: Rq, R^^ = X^, Xg (EX) Rg, Rg PRESERVED 

EXECUTION TIME: 13R + W + 51N 



TEMPORARY LOCATIONS USED: N OOEO 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 1 FOR INPUT EXPONENT ONLY 
(MAXIMUM DEPTH) 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

2. WARNING: ACCESS THIS SUBROUTINE ONLY BY A JSR INSTRUCTION! 
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SUBROUTINE: LEFT NORMALIZE 



LABEL: LFNOR 



DESCRIPTION: LEFT NORMALIZES A FLOATING-POINT NUMBER IN R^, R;^ 



ENTRY POINT: FC97 

ENTRY CONDITIONS: Rq, R^ = X 

EXIT CONDITIONS: Rq, Rj^ = X NORMALIZED 

EXECUTION TIME: 449R + 24W + 1745N 



;R2, R3 ARE PRESERVED 



TEMPORARY LOCATIONS USED: OOEO, OOEl 

SUBROUTINES CALLED: CZERO, FPCOMP, EXTEXP, DSHL, ADDEXP 

STACK WORDS NEEDED: 2 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: DOUBLE LEFT NORMALIZE 



LABEL: DLNORM 



DESCRIPTION: NORMALIZES TWO FLOATING-POINT NUMBERS Rq, R-^ AND Rg, Rg, WITH THE 

RESULTANTS SWITCHED (IN REGISTERS Rg, Rg AND Rq, R^^, RESPECTIVELY). NEW 
EXPONENTS ARE IN MEMORY. 



ENTRY POINT: FC8A 

ENTRY CONDITIONS: Rq, Rj^ = X^, Xg (EX) Rg, Rg = Y^^, Y^ (EY) 

EXIT CONDITIONS: Rq, Rj^ = Y NORMALIZED (EX = 0) Rg, R3 = X NORMALIZED (EX = 0) 

00E6 = EY 00E5 = EX 

EXECUTION TIME: 947R + 52W + 3691N 



TEMPORARY LOCATIONS USED: 00E0-00E2, 00E5, 00E6 

SUBROUTINES CALLED: LFNOR, EXTEXP 

STACK WORDS NEEDED: 3 
(MAXIMUM DEPTH) 



COMMENTS: 



ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: FRACTIONAL TO FLOATING POINT 



LABEL: FLOAT 



DESCRIPTION: CONVERTS DOUBLE -PRECISION FRACTIONAL NUMBER TO FLOATING-POINT 
FORMAT BY TRUNCATING LEAST SIGNIFICANT 8 BITS. 



ENTRY POINT: FC85 

ENTRY CONDITIONS: R^, R^^ = X-^, X2 

EXIT CONDITIONS: R^, R^^ = X^, Xg' (0) 



EXECUTION TIME: 



3R + 9N 



Rg, R3 PRESERVED 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: NONE 
(MAXIMUM DEPTH) 



COMMENTS: 



ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: FLOATING POINT TO FRACTIONAL 



LABEL: SFO 



DESCRIPTION: CONVERTS A FLOATING-POINT NUMBER TO A FRACTIONAL DOUBLE -PRECISION 
NUMBER 



Rf 



X, 



MANTISSA 



Ri 



^2 


EX 



EXPONENT 



BEFORE 



X 



X 



DOUBLE PRECISION 



AFTER 



ENTRY POINT: FCC 7 

ENTRY CONDITIONS: Rq, R-^ = X^, Xg (EX) 

EXIT CONDITIONS: Rg, R^ = X DOUBLE PREC. 

EXECUTION TIME: 207R + W + 822N 



R2, R3 PRESERVED 



TEMPORARY LOCATIONS USED: OOEO 

SUBROUTINES CALLED: EXTEXP, DSH 

STACK WORDS NEEDED: 2 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

CAUTION: VALUE OF INPUT MUST BE < 1. IN MAGNITUDE. 
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National Semiconductor Corporation 

2900 Semiconductor Drive 
Santa Clara, California 95051 
(408) 732-5000 
TWX: 910-339-9240 



National Semiconductor GmbH 

D 808 Fuerstenfeldbruck 
Industrrestrasse 10 
West Germany 
Telephone: (08141) 1371 
Telex: 05-27649 



NS Electronics SDN BHD 

Batu Berendam 

Free Trade Zone 

Malacca, Malaysia 

Telephone: 5171 

Telex: NSELECT 519 MALACCA (c/o Kuala Lumpur) 



National Semiconductor (UK) Ltd. 

Larkfield Industrial Estate 
Greenock, Scotland 
Telephone: GOUROCK 33251 
Telex: 778 632 



NS Electronics (PIE) Ltd. 

No. 1100 Lower Delta Rd. 
Singapore 3 
Telephone: 630011 
Telex: NATSEMI RS 21402 



REGIONAL AND DISTRICT SALES OFFICES 



ALABAMA 

DIXIE REGIONAL OFFICE 
3322 Memorial Parkway, S.W. #67 
Huntsville, Alabama 35802 
(205) 881-0622 
TWX: 810-726-2207 

ARIZONA 

*ROCKY MOUNTAIN REGIONAL OFFICE 
7349 Sixth Avenue 
Scotlsdale, Arizona 85251 
(602) 945-8473 
TWX: 910-950-1195 

CALIFORNIA 

*NORTH-WEST REGIONAL OFFICE 
2680 Bayshore Frontage Road, Suite 112 

/Mountain View, California 94043 
(415) 961-4740 
TWX: 910-379-6432 

NATIONAL^EMICONDUCTOR 
*LOS ANGELES REGIONAL OFFICE 

Valley Freeway Center Building 

15300 Ventura Boulevard, Suite 305 

Sherman Oaks, California 91403 

(213) 783-8272 

TWX: 910-495-1773 

NATIONAL SEMICONDUCTOR 
*SOUTHERN CALIFORNIA REGIONAL OFFICE 

17452 Irvine Boulevard, Suite M 

Tustin, California 92680 

(714) 832-8113 

TWX: 910-595-1523 

CONNECTICUT 

DISTRICT OFFICE 
Commerce Park 
Danbury, Connecticut 06810 
(203) 744-2350 
'MID-ATLANTIC REGIONAL SALES OFFICE 
25 Sylvan Road South 
Westport, Connecticut 06880 
(203) 226-6833 



FLORIDA 

'CARIBBEAN REGIONAL OFFICE 
2721 South Bayshore Drive, Suite 121 
Miami, Florida 33133 
(305) 446-8309 
TWX: 810-848-9725 

ILLINOIS 

NATIONAL SEMICONDUCTOR 
WEST-CENTRAL REGIONAL OFFICE 
800 E. Northwest Highway, Suite 203 
Mt. Prospect, Illinois 60056 

(312) 394-8040 
TWX: 910-689-3346 

INDIANA 

NATIONAL SEMICONDUCTOR 
NORTH-CENTRAL REGIONAL OFFICE 
P.O. Box 40073 
Indianapolis, Indiana 46240 
(317)255-5822 

MARYLAND 

CAPITAL REGIONAL SALES OFFICE 
300 Hospital Drive, No. 232 
Glen Burnie, Maryland 21061 
(301) 760-5220 
TWX: 710-861-0519 

MASSACHUSETTS 

*NORTH-EAST REGIONAL OFFICE 
#8WallisCt. 

Lexington, Massachusetts 02173 
Telephone: (617) 861-6090 
TWX: 710-332-0166 

MICHIGAN 

'REGIONAL OFFICE 
23629 Liberty Street 
Farmington, Michigan 48024 

(313) 477-0400 
TWX: 810-242-2902 



MINNESOTA 

DISTRICT SALES OFFICE 

8053 Bloomington Freeway, Suite 101 

Minneapolis, Minnesota 55420 

(612)888-3060 

Telex: 290-766 

NEW JERSEY 

DISTRICT SALES OFFICE 

140 Sylvan Avenue 

Englewood Cliffs, New Jersey 07632 

(201) 461-5959 

TWX: 710-991-9734 

AREA OFFICE 

14 Commerce Drive 

Cranford, New Jersey 07016 

(201 ) 272-3344 

TWX: 710-996-5803 

NEW YORK 

CAN-AM REGIONAL SALES OFFICE 
104 Pickard Drive 
Syracuse, New York 13211 
(315) 455-5858 

OHIO 

DISTRICT SALES OFFICE 
Financial South Building 
5335 Far Hills, Suite 214 
Dayton, Ohio 45429 
(513)434-0097 
TWX: 810-459-1615 

TEXAS 

'SOUTH-CENTRAL REGIONAL OFFICE 
5925 Forest Lane, Suite 205 
Dallas, Texas 75230 
(214) 233-6801 
TWX: 910-860-5091 

WASHINGTON 

DISTRICT OFFICE 
300 120th Avenue N.E. 
Building 2, Suite 205 
Bellevue, Washington 98005 
(206) 454-4600 



INTERNATIONAL SALES OFFICES 

AUSTRALIA 

NS ELECTRONICS PTY, LTD. 

Cnr. Stud Road & Mountain Highway 

Bayswater, Victoria 3153 

Australia 

Telephone: 03-729-6333 

Telex: 32096 

BELGIUM 

NATIONAL SEMICONDUCTOR BELGIUM 
789 Ave. Houba de Strooper 
1020 Bruxelles 
Telephone: 02-478-3400 
Telex: 61 007 NatSem B 

CANADA 

NATIONAL SEMICONDUCTOR 
DISTRICT OFFICE 
1111 Finch Avenue West 
Suite 154 

Downview, Ontario 
M3J 2E5 
(416) 630-5751 
TWX: 610-492-1337 

DENMARK 

NATIONAL SEMICONDUCTOR DENMARK 

Nyhavn 69 

1051 Copenhagen 

Telephone: (1) 153110 

Telex: 160 39 

*Wicroproces8or System Specialist Available 



ENGLAND 

NATIONAL SEMICONDUCTOR (UK) LTD. 

The Precinct 

Broxbourne, Hertfordshire 

EN 107 HY 

England 

Telephone: Hoddesdon 69571 

Telex: 267-204 

FRANCE 

NATIONAL SEMICONDUCTOR FRANCE 
EXPANSION 10000 
28 rue de la Redoute 
92-260 Fountenay Aux Roses 
Telephone: 660.81.40 
Telex: NSF 25956F-H 

GERMANY 

NATIONAL SEMICONDUCTOR GmbH 
8000 Munchen 81 
Cosimastr. 4/1 
Telephone: 089/915027 
Telex: 05-22772 

HONG KONG 

NS ELECTRONICS (HONG KONG) Ltd. 

11th Floor 

4 Ming Yip Street 

Kwun Tong 

Kowloon, Hong Kong 

Telephone: 3-411241-8 

Telex: 73866 NSE HK HX 

Cable: NATSEMI 



ITALY 

NATIONAL SEMICONDUCTOR ITALY 
Via Valassina 24 , 
20159 Milano 
Telephone: (02)688 4617 
Telex: 36-540 

JAPAN 

'NATIONAL SEMICONDUCTOR JAPAN 
Nakazawa Building 
1-19 Yotsuya, Shinjuku-Ku 160 
Tokyo, Japan 
Telephone: 03-359-4571 
Telex: J 28592 

SWEDEN 

NATIONAL SEMICONDUCTOR SWEDEN 

Sikvagen17 

13500 Tyreso-Stockholm 

Telephone: 08/7 1204 80 

Telex: 11293 

TAIWAN 

NS ELECTRONICS (HK) LTD. 
TAIWAN LIAISON OFFICE 
#60 Teh Hwei Street 
P.O. Box 68-332 
Taipei Taiwan ROC 
Telephone: 553354 
Cable; NSTW TAIPEI 
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